Modeling elastic and photoassisted transport in organic molecular wires: length 
dependence and current-voltage characteristics 
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Using a 7r-orbital tight-binding model, we study the elastic and photoassisted transport prop- 
erties of metal-molecule-metal junctions based on oligophenylenes of varying lengths. The effect 
of monochromatic light is modeled with an ac voltage over the contact. We first show how the 
low-bias transmission function can be obtained analytically, using methods previously employed for 
simpler chain models. In particular, the decay coefficient of the off-resonant transmission is ex- 
tracted by considering both a finite-length chain and infinitely extended polyphenylene. Based on 
these analytical results, we discuss the length dependence of the linear-response conductance, the 
thermopower, and the light-induced enhancement of the conductance in the limit of weak intensity 
and low frequency. In general the conductance-enhancement is calculated numerically as a function 
of the light frequency. Finally, we compute the current-voltage characteristics at finite dc voltages, 
and show that in the low-voltage regime, the effect of low-frequency light is to induce current steps 
with a voltage separation determined by twice the frequency. These effects are more pronounced 
for longer molecules. We study two different profiles for the dc and ac voltages, and it is found 
that the results are robust with respect to such variations. Although we concentrate here on the 
specific model of oligophenylenes, the results should be qualitatively similar for many other organic 
molecules with a large enough electronic gap. 

PACS numbers; 73.50.Pz,85.65.+h,73.63.Rt 



o 

> 

m 
m 



o 

00 

o 



X 



I. INTRODUCTION 



The use of single-molecule electrical contacts for op- 
toelectronic purposes such as light sources, light sen- 
sors, and photovoltaic devices is an exciting idea. Yet, 
due to the difhculties that light-matter interactions 
in nanoscale systems pose for theoretical and exper- 
imental investigations, the possibilities remain largely 
unexplored. Concerning experiments, it has been 
shown that Hght can be used to change the conforma- 
tion of some molecules even when they are contacted 
to metallic electrodes, thus enabling light-controlled 
switching.^ Some evidence of photoassisted processes in- 
fluencing the conductance of laser-irradiated metallic 
atomic contacts has also been obtained.^ Theoretical in- 
vestigations of light-related effects in molecular contacts 
are more numerous^^^^^i^^^iii^i^^i^^ii^ia but 
they are mostly based on highly simplified models, 
whose validity remains to be checked by more de- 
tailed calculations^^'^^ and experiments. However, for 
the description of the basic phenomenology, model ap- 
proaches can be very fruitful, as they have been in stud- 
ies of elastic transport in the past. Properties of lin- 
ear single-orbital tight-binding (TB) chains, in partic- 
ular, have been studied in detail, and to a large part 
analytically^^iS3^'25^i27^^^^^ In a step towards 
a more realistic description of the geometry, symmetries, 
and the electronic structure of particular molecules, em- 
pirical TB approaches such as the (extended) Hiickel 



method have proved useful. 



4.8.33.34.35 



Based on a combination of density-functional calcu- 
lations and simple phenomcnological considerations, we 
have recently described the photoconductance of metal- 
oligophenylene-metal junctions!^ It was discussed how 
the linear-response conductance may increase by orders 
of magnitude in the presence of light. This effect can 
be seen as the result of a change in the character of the 
transport from off-resonant to resonant, due to the pres- 
ence of photoassisted processesi^i^ Consequently, the 
decay of the conductance with molecular length is slowed 
down, possibly even making the conductance length- 
independent i^i^ 

In this paper we apply a Hiickel-type TB model of 
oligophenylene-based contacts^^ combined with Green- 
function methods^ to study the effects of monochromatic 
light on the dc current in metal-oligophenylene-metal 
contacts. Again we concentrate on the dependence of 
these effects on the length of the molecule. We begin 
with a detailed account of the elastic transport prop- 
erties of the model, and show that the zero-bias trans- 
mission function can be obtained analytically, similarly 
to simpler chain models;—!^ We demonstrate how infor- 
mation about the length dependence of the transmission 
function for a finite wire can be extracted from an in- 
finitely extended polymer. Based on these analytical re- 
sults, we discuss the length dependences of the conduc- 
tance and the photoconductance for low-intensity and 
low- frequency light. While the conductance decays ex- 
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ponentially with length, its relative enhancement due to 
light exhibits a quadratic behavior. Here we also briefly 
consider the thermopower, whose length dependence is 
linear. Next, we calculate numerically the zero-bias pho- 
toconductance as a function of the light frequency uj, 
and find that the conductance-enhancement due to light 
is typically very large. In particular, we show that 
the results of Ref. 5 are expected to be robust with re- 
spect to variations in the assumed voltage profiles. Fi- 
nally, we describe how the step-like current-voltage (/- 
V) characteristics are modified by light. At high lo the 
most obvious effect is the overall increase in the low- 
bias current. At low w, additional current steps similar 
to those in microwave-irradiated superconducting tunnel 
junctions*^^'^^ can be seen. Their separation, in our case 
of symmetric junctions, is roughly 2TiLo/e. 

TB models of the type we shall consider neglect vari- 
ous interaction effects (see Sec. |V] for a discussion), and 
thus cannot be expected to give quantitative predictions. 
However, the qualitative features of the results rely only 
on the tunneling-barrier character of the molecular con- 
tacts, which results from the fact that the Fermi energy 
of the metal lies in the gap between the highest-occupied 
and lowest-unoccupied molecular orbitals (HOMO and 
LUMO) of the molecule. Thus, these features should re- 
main similar for junctions based on many other organic 
molecules exhibiting large HOMO-LUMO gaps. The 
light-induced effects, if verified experimentally, could be 
used for detecting light, or as an optical gate (or "third 
terminal") for purposes of switching. 

The rest of the paper is organized as follows. In Sec.HIl 
we describe our theoretical approach, discuss the gen- 
eral properties of TB wire models, and introduce the 
Green-function method for the calculation of the elas- 
tic transmission function. Then, in Sec. IIIII we calculate 
the transmission function of oligophenylene wires analyt- 
ically. The decay coefficient for the off-resonant transmis- 
sion is extracted also from infinitely extended polypheny- 
lene. Following that, in Sec. II VI we present our numerical 
results for the conductance, the thermopower, the pho- 
toconductance, and the I-V characteristics. Finally, Sec. 
IVl ends with our conclusions and some discussion. De- 
tails on the calculation of the time-averaged current in 
the presence of light are deferred to the appendixes. In 
App.|A]a simplified interpretation of the current formula 
is derived, and in App. |B] a brief account of the general 
method is given. Readers mainly interested in the discus- 
sion of the results for the physical observables can skip 
most of Sees. HI] and Uni and proceed to Sec. IIVI 



II. THEORETICAL FRAMEWORK 

A. Transport formalism 

Our treatment of the transport characteristics for the 
two-terminal molecular wires is based on Green's func- 
tions and the Landauer-Biittiker formalism, or its gen- 



eralizations. Assuming the transport to be fully elastic, 
the dc electrical current through a molecular wire can be 
described with 

I{V) = ^ 1 dET{E, V)[fL{E) - ME)]. (1) 

Here V is the dc voltage and t{E, V) is the voltage- 
dependent transmission function, while fx{E) = 
l/[exp((£^ — fJ-x)/kBTx) + 1], fJ-x, and Tx are the Fermi 
function, the electrochemical potential, and the temper- 
ature of side X = L,R, respectively.^^ The electrochem- 
ical potentials satisfy eV — Afi — fj.^ — fj,ji, and we 
can choose them symmetrically as fij^ = Ep -j- eV/2 and 
/i/j = Ep — eV/2, where Ep is the Fermi energy. For 
studies of dc current we always assume = — 0. 
Of particular experimental interest is the linear-response 
conductance Gdc — dl /dV\v=o, given by the Lan- 
dauer formula Gdc — Got{Ep), where Go = 2e^//i and 
t{E) — t{E, V — 0). In most junctions based on organic 
oligomers, the transport can be described as off-resonant 
tunneling. This results in the well-known exponential 
decay of Gdc with the number N of monomeric units in 
the molecule4S At finite voltages V, the current increases 
in a stepwise manner as molecular levels begin to enter 
the bias window between fj.^ and fifi (Ref. |24|). We shall 
consider both of these phenomena below. 

If a small temperature difference AT — Tl — Tn at 
an average temperature T = [T^ + Tfi)/2 is applied, 
heat currents and thermoelectric effects can arise i^^'^^i"^^ 
In an open-circuit situation, where the net current / 
must vanish, a thermoelectric voltage A/z/e is generated 
to balance the thermal diffusion of charge carriers. In 
the linear-response regime the proportionality constant 
S — — (A/i/eAT)/=o is the Seebeck coefficient. We will 
briefly consider this quantity below as an example of an 
observable with a linear dependence on the molecular 
length N , but will not enter a more detailed discussion 
of thermoelectricity or heat transport. 

The quantity we are most interested in is the dc cur- 
rent in the presence of monochromatic electromagnetic 
radiation, which we refer to as light independently of its 
source or frequency uj. We model the light as an ac volt- 
age with harmonic time-dependence V{t) — VacCOs{ujt) 
over the contact. The current averaged over one period 
of V{t) can be written in the form^iii^ 

I{V-a,u;)=- J2 dE[T^^l{E,V-a,u) 

k=-oo-^ (2) 

X hiE) ~ Ti''liE,V;a,Lo)fR{E)]. 

Here the transmission coefficient T^j^{E), for example, 
describes photoassisted processes taking an electron from 
left [L) to right (i?), under the absorption of a total of 
k photons with energy Tiuj. The parameter a = eVac/TT-uj 
describes the strength of the ac drive4i It is deter- 
mined by the intensity of the incident light and possi- 
ble field-enhancement effects taking place in the metallic 
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FIG. 1: (Color online) A finite block chain of length N=3 
connected to electrodes at its two ends. This gives rise to self 
energies Sn and Sjvjv on the terminating blocks. 



nanocontactJ^ Again, in addition to the full I-V char- 
acteristics, we study in more detail the case of linear 
response with respect to the dc bias, i.e., the photocon- 
ductance Gdc{Q:,oj) = dI{V; a,LLj)/dV\v=o- The argu- 
ments a and uj distinguish it from the conductance Gdc, 
although we sometimes omit a for notational simplicity. 

The calculation of the coefficients T^l/Lf({E) is rather 

complicated in general)^ and we defer comments on this 
procedure to App. [Bl Below we shall mostly refer to an 
approximate formula (see App. [X| that can be expressed 
in terms of t{E). This amounts to a treatment of the 
problem on the level of the Tien-Gordon approach i^i^i^ 
The full Green-function formalism for systems involving 
ac driving is presented in Ref. |3- 

In noninteracting (non-self-consistent) models it is in 
general not clear how the voltage drop should be di- 
vided between the different regions of the wire, and 
the electrode-wire interfaces. A self-consistent treatment 
would be in order in particular for asymmetrically cou- 
pled molecules. We only concentrate on left-right sym- 
metric junctions, where where both the dc and ac volt- 
ages {V and Vac) are assumed to drop according to one 
of two different symmetrical profiles. The symmetry of 
the junctions excludes rectification effects, such as light- 
induced dc photocurrents in the absence of a dc bias 
voltage . ^'^'^^ However, light can still have a strong in- 
fluence on the transmission properties of the molecular 
contact, as will be discussed below. It will be shown 
that our conclusions are essentially independent of the 
assumed voltage profile. 



B. Wire models 

Below we will specialize to the case of a metal- 
oligophenylene-metal junction. However, to make some 
general remarks, let us first consider a larger class of 
molecular wires that can be described as N separate 
units forming a chain, where only the nearest neighbors 
are coupled (see Fig. [J). We only discuss the calcula- 
tion of the elastic transmission function t{E, V) here, as 
this will be the focus of our analytical considerations in 
Sec. mil From this quantity (at V — 0), the various 
linear-response coefficients such as the conductance and 
the thermopower can be extracted. Furthermore, as al- 
ready mentioned, it suffices for an approximate treatment 
of the amplitudes t^j^I{E) as well. 



We assume a basis \xp^^) of local (atomic) orbitals, 
where p = 1, . . . , N indexes the unit, while a = 1, . . . , Mp 
denotes the orbitals in each unit4i For simplicity, the 
basis is taken to be orthonormal, i.e. iXp^'^lxf') = 



^afs^pq- The (time- independent) Hamiltonian Hj 



(a:/3) 



-Jpq 



form 



of the wire is then of the block-tridiagonal 



/Hii Hi2 
H21 H22 



H 



H 



23 



Hn-1,N-2 -ffjV-l,JV- 



1 -f^TV-l.JV 

Hnn / 



(3) 



where Hpg with p,q = 1, . . . ,N are Mp x Mq matrices. 
(The unindicated matrix elements are all zeros.) 

In the non-equilibrium Green-function picture, the ef- 
fect of coupling the chain to the electrodes is described in 
terms of "lead self energies" We assume these to be lo- 
cated only on the terminal blocks of the chain, with com- 
ponents Sii and 'Snn- The inverse of the stationary- 
state retarded propagator for the coupled chain will then 
be of the form 



/-Fii hi2 

h21 h22 



h 



\ 



23 



h 



V 



Ar-l,Af-2 



Fnn / 



(4) 



H 



hr, 



\pp - 



- Hpp , and 



Here hp,p±i 

E+iO"^, while Fu — /in — En and Fnn = hNN — '^NN- 
Charge-transfer effects between the molecule and the 
metallic electrodes shift the molecular levels with respect 
to the Fermi energy Ep- In a TB model, these can be rep- 
resented by shifting the diagonal elements of H. Once a 
transport voltage V is applied, further shifts are induced. 
In our model the voltage-induced shifts will be taken from 
simple model profiles, and the relative position of Ep will 
be treated as a free parameter. 

Effective numerical ways of calculating the propagator 
G = F^^ for block-tridiagonal Hamiltonians exist 
In Sec, mil we shall be interested in a special case, where 
Hp^p-i — H i, Hp,p+i = Hi and Hpp = Hq with 
the same Hi ~ H'^-^ and Hq (of dimension Mp ~ M) 
for all p, describing an oligomer of identical monomeric 
units. In such cases also analytical progress in calculating 
the current in Eq. ([T]) may be possible. Once the Green 
function G is known, the transmission function is given 
by48 



r(S, V) = Tr[rnGiJvrArjv(GiAr)t], 



(5) 



where Tn = -2ImSii and T,ii{E,V) = Y,ii{E - eV/2), 
for example. 

Typically Ep lies within the HOMO-LUMO gap, re- 
sulting in the exponential decay t[Ep) ~ ^-Pi^p)^ with 



TV, characteristic of ofF-resonant transport. The decay co- 
efficient P{Ep) is actually independent of Sn and Sat at. 
This can be seen by considering the Dyson equation 
G = Q + ^SG, where G and Q are the Green function 
of the coupled and uncoupled wires, respectively, and S 
is the matrix for the lead self-energies. Assuming that 
QiN decays exponentially with N, then 



,(a) 



,(a) 



,(a) 



GiAT sa (1 - giiSii) ^Qi 



N 



(6) 



when N —^ oo, and therefore Gin decays with the same 
exponent. Thus, one can in principle obtain the decay 
exponent from the propagator of an isolated molecule, or 
even an infinitely extended polymer. In the next Section 
we demonstrate this by extracting the decay exponent 
of a finite oligophenylene junction from the propagator 
for polypheny lene. We note that in doing so, we neglect 
the practical difficulty of determining the correct relative 
position of Ep. 

There are efficient numerical methods for computing 
the lead self-energies for different types of electrodes and 
various bonding situations between them and the wire. 
Typically, the methods are based on the calculation of 
surface Green's functionsi^ Below we shall simply treat 
the self-energies as parameters. 



III. PHENYL-RING-BASED WIRES 

In this Section we discuss a special case of the type 
of wire model introduced above, describing an oligomer 
of phenyl rings coupled to each other via the para (p) 
position:^ The bias voltage V is assumed to be zero. 
In the special case that we will consider, the inver- 
sion of Eq. ^ can then be done analytically with the 
subdeterminant method familiar from elementary linear 
algebraj22i^i22i^ Below, we first use this method for cal- 
culating the propagator of the finite-wire junction and 
derive the decay exponent (3{E) of the transmission func- 
tion at off-resonant energies. After that we rederive 
the decay exponent by considering an infinitely extended 
polymer of phenyl rings. 




FIG. 2: (Color online) A finite chain of length ^^=3 connected 
to electrodes at its two ends. This gives rise to self energies 
Sii and Sivjv on the end sites. The nearest-neighbor hop- 
pings inside the ring (—7) and between the rings (—17) are 
different. The lower part indicates also the numbering of the 
M — 6 carbon atoms within a ring. 



The ring-tilt angle (p can be controlled to some ex- 
tent using side groups. For example, two side groups 
bonded to adjacent phenyl rings can repel each other ster- 
ically, thus increasing the corresponding tilt anglci^i^ In 
fact, even the pure oligophenylenes in the uncharged state 
have (fi — 30° — 40° due to the repulsion of the hydro- 
gen atoms. However, side groups can introduce also 
"charging" or "doping" effects, which shift the molecular 
levels 1^ 

For definiteness, we number the M = 6 carbon atoms 
of a phenyl ring according to the lower part of Fig. O The 
corresponding orbitals appear in the basis in this order. 
Thus the blocks in Eq. ([3]) are 
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-7 



9:9 






V 

ior q ^ 1, . . . , N and 



-7 

,(2) 
^9 



"7 



J3) 



-7 







.(5) 
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7 -7 



o\ 




-7 

-7 
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(7) 



A. Oligo-p-phenylene junction 

Our model for the oligophenylene-based molecular 
junction is depicted in Fig.[2l Within a simple 7r-electron 
picture, the electronic structure of the oligophenylene 
molecule can be described with a nearest-neighbor TB 
model with two different hopping elements —7 and —rj 
(Ref. 52). Here —7 is for hopping within a phenyl ring, 
between the p orbitals oriented perpendicular to the 
ring plane, while —77 describes hopping between adjacent 
rings. Due to the symmetry of the orbitals, the mag- 
nitude of 77 depends on the ang le if between the rings 
proportionally to cos(^ (Ref. [531). We shall assume that 
r] = 7C0S(/3, and thus I??] < 7. In this way the natural 
energy scale of the model is set by 7 alone. 
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(8) 



with Hq^i_q ~ [Hq^g+i]'^. Hcrc the onsite energies €q 
may be shifted non-uniformly to describe effects of pos- 
sible side-groupsi^ For simplicity, we shall consider all 
phenyl rings to have a similar chemical environment, and 
thus all onsite energies are taken to be equal. 

As a first step we note that, assuming e^"^ = eq for 
all a, the eigenvalues for the Hamiltonian Hqq of the 
isolated unit are e,— 7, eq-l-7, e,— 7, eq+j, e, — 27, eq-l-27. 
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while the corresponding orthonormahzed eigenvectors are 
-^(0, -1, 1, -1, 1, of, -i=(0, 1, -1, -1, 1, 0)^, 

'=(-2,-l,-l,M,2f,^(2,-l,-l,-l,-l,2f, 



(9) 



-^(1, 1, 1, 1, 1, 1)^, -i=(-l, 1, 1, -1, -1, 1)^. 



The first two of the eigenstates have zero weight on the 
ring-connecting carbon atoms 1 and 6. Therefore, these 
eigenstates do not hybridize with the levels of the ad- 
jacent rings and consequently cannot take part in the 
transport. This will be seen explicitly in the derivation 
of the propagator. We note that these results can also 
be used to determine a realistic value for the hopping 7 
from the HOMO-LUMO splitting of benzene^ 

Below we shall only consider the analytically solvable 
case, where all onsite energies are set to the same value. 
We choose this value as our zero of energy: e,"^ = for all 
q = 1, . . . ,N and a = 1, . . . , M. Later on we shall relax 
this assumption in order to describe externally applied dc 
and ac voltage profiles. In the absence of such voltages, 
the inverse propagator [Eq. (g])] consists of the blocks 
hp^p — ho, hp^p^i — h^i, and /ip,p+i = hi, where 
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(10) 



h , = 



/O 77\ 









Vo oy 



and hi = [h^i]^ . The leads are assumed to couple 
only to the terminal carbon atoms, thus making the self- 
energies 6x6 matrices of the form 












/O \ 



NN 
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\ooooy yo---o Y.r) 

(11) 

We also define the symbol "tilde" ("'), which means the 
replacement of the first column of a matrix by 77 followed 
by zeros. For example 
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(12) 



For the evaluation of Eq. ([5]), we only need the com- 
ponent Gi^MN = [Gin]im- Using the subdeterminants 
of F G^^, we have 



Gi 



MN 



(-l)^^-^^+idet[F(MiV|l)] 
det[F] ■ 



(13) 



Here 0{i, . . . , k\j, . . . ,1) is the submatrix of O obtained 
by removing the rows i, . . . ,k and columns j, . . . ,1. We 
shall also denote by L and R the "leftmost" and "right- 
most" row or column of a matrix. Thus, for example 
det[F{MN\l)] = det[F{R\L)] 

Let us first concentrate on the denominator of Eq. 
(fT5)) . It is easy to see that det[F] can be written in 
terms of determinants related to the inverse Green func- 
tion ^ — Q^^ of the uncoupled wire as foUows'^'^ 

det[F] = dct[:F] - Sl det[:F(i|i)] - T.Rdct[J^{R\R)] 
+ SL^Rdet[J^{L,R\L,R)]. 

(14) 

Furthermore, due to the symmetry of the molecule, 
det[J^{R\R)] = det[:F(L|L)]. Thus we are left with calcu- 
lating three types of determinants. It can be shown that, 
for 1 < n < iV, all of them satisfy a recursion relation of 
the form 



= iEl-^^)Y 



{El 



1') 



£)("-!) 

a -c\ 
c b 



(15) 



For example, in the calculation of det[^], we have D'^") = 

det[:F("'] and = det[i^^"''], where the additional 

superscript (n) on the matrices denotes the number of 
the M X M diagonal blocks. The elements of the matrix 
Y are given by 



h=-rf{El-^^) 
c = vE+iEl-3j^) 



■47^) 



(16) 



Only the initial condition (?i = 1) and the last step of the 
recursion (n — N) will differ for the three determinants. 
The recursion relations can be solved by calculating Y" 
explicitly, which can be done by diagonalizing Y. The 
eigenvalues of Y are Ai,2 = {a + b^ ^^(a — by — 4c2)/2, 
while the (unnormalized) eigenvectors are 



Vl,2 = 



bT y/{a - by - Ac^ 



(17) 



Then, if V = {vi,V2) and A 
yn ^ vM'V-^. The result is 



diag(Ai,A2), we have 



Y'" 



(«) (n) 



(18) 



where the components are given by 



(«) _ (A? " Ag)(b ~ fl) + {VI + A^')v/(a - 6)2 ~ 4^^ 
2^ (a - 6)2 - 4c2 



(„) (A? - A^)(a - 6) + (A';' + \'i)^{a - 6)2 - 4c2 

2/22 = 



2^(a - 6)2 - 4c2 



(n) _ in) _ C(A^ - A^) 

y\i — 2/21 



^(a - 6)2 - 4c2 



(19) 



Using these, we can now write exphcit expressions for the 
three required determinants. For det[.^], the recursion 
can be started at n = 1 with the initial conditions Z?*^"^ = 
1 and 1)^°-' = and carried out up to n = N . The resuh 
is 



det[^W] = (i?2 _^2yv^.» 



(20) 



The other two determinants require special initial and 
final steps, and the results are 



det[^W(L|L)]=(i?^-7^)^^fV^ 
det[JPW(L,i?|i,i?)]=(i?2 _^2)^[y(^-i)^ 



(21) 



Next, we consider the determinant in the numerator 
of Eq. (Hll), det[F(^'(i?|L)] = Aq\\T^^'\R\L)\. It can 
easily be shown that it satisfies the recursion relation 

det[:rW(i?|L)] = 2^-{'[EX--1'')A(ti{T^^-^\R\L)\ 

(22) 

and so 

dct[^W(i?|i)] = 2^(^7')^(£;2 _ ^2^1^, (23) 
Now, the Green function of Eq. IT^l) can be written as 



Gi.MN — 



2/11 + ^Li?2/21 /'? + ^L^-R. (^2/21 C - 7/22 Oj 

(24) 

where we used the shorthand Sla — + S_r. 

It is notable that the common (iJ^ ~7^)^ factors can- 
celed out from the final propagator. These factors ap- 
parently correspond to the two eigenvectors of Hq [Eq. 
([5])] having zero weight on the ring-connecting atoms 1 
and 6. The cancellation is a manifestation of the physical 
fact that such localized states cannot contribute to the 
transport through the molecule. In the infinite polymer 
to be discussed below, these states appear as completely 
flat bands in the band structure. 

To conclude this part, we point out that for E inside 
the HOMO-LUMO gap [more precisely, when (a - 6)2 - 
4c2 > 0] the eigenvalues A 1,2 are real- valued and the de- 
cay exponent of the transmission t{E) for large N is con- 
trolled by the one with a larger absolute value. Since 



(a) 




(b) 
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FIG. 3: (Color online) Phenyl-ring chains: (a) a periodic 
chain with A'^ units and (b) an infinite chain. Case (b) is 
obtained from (a) in the limit A'^ 00. 



inside the gap £' sa 0, we find that A2 > Ai > 0. Then, 
using Eq. ([5]) and omitting iV- independent prefactors, the 
decay of the transmission for large N follows the law 



-{E) 



\2{E) 



-2N 



_ -2N\n[\2(E)/{2iri'')] 



Thus the decay exponent is given by 

/3(i?) = 21n[A2(ii;)/(2ry73)]. 



(25) 



(26) 



We note that for resonant energies, oscillatory depen- 
dence of t{E) on N can be expected, instead, and for 
limiting cases also power-law decay is possibles^ Next, 
we shall reproduce the result for the decay exponent by 
considering an infinitely extended polymer. 



B. Poly-p-phenylene 

For comparison with the "correct" evaluation of the 
propagator and the decay coefficient for a finite chain, 
let us consider the propagator for an infinitely extended 
polymer. To describe the polymer, we start from a finite 
chain with periodic boundary conditions. Neglecting cur- 
vature effects, the latter actually represents a ring-shaped 
oligomer, as depicted in Fig.[3]Ja). 

Let us first consider the eigenstates of the periodic 



chain. The Hamiltonian H, 



(",/3) 



the general form 



pq 



{xt^mx'f') is of 



H 



Hi Ho Hi 
\Hi H i Ho J 



(27) 



where ifo,±i are the M x M matrices (M = 6) of Eqs. 
O and dHD, with e^"^ = 0. (Again, only nonzero elements 
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are indicated.) The normalized eigenvectors ■4>^^\k) sat- 
isfying 



(28) 



are of the Bloch form -0^"' (fc) = e*'^'«'' (fc) /a/TV, where 
arc the normaUzed eigenvectors of 

H{k) = e'^'^Hi + Ho + e-'^'^H^i (29) 

with the eigenvalue E^"\k), and n — 1, . . . , M. Due to 
the finiteness of the wire, the k values are restricted to 
kfj, = 2'Kfj,/Nd, where /i is an integer and d is the lattice 
constant (the length of a single phenyl-ring unit). 

The spectral decomposition of the (retarded) propaga- 
tor g{E) = (i?+l — H)^^ of the chain is of the form 



9, 



{E) 



E^ 



E(^){k^ 



(30) 



with the Bloch states 
1 



L7V/2J 

N ^ 

p=-rAr/2l + l 



M 
a=l 



(31) 

In the limit of large N [Fig. W}^)] , we can use N ^ J2^i ~^ 
{d/2'K) /^^^^ dfc to turn the summation into an integral 
over the first Brillouin zone. In this case, there are M = 6 
bands with energies 



±^V'f + 572-2i?(fc) (32) 
1 



v2 



where 



Clearly we have the symmetries E^^''{k) = — i?^^)(fc), 
i;(3)(A;) = -£;(4)(fc), and i;(5)(fc) = -E'-^^k). For n = 
1, 2 the bands are completely flat, and the corresponding 
eigenvectors <p^^''^\k) are as in Eq. i.e., independent 
of k and completely localized on atoms a = 2,3,4,5. 
Thus for p q, they do not contribute to the propaga- 
tor in Eq. ((30| . For n = 3,4,5,6, the vectors are very 
complicated, but they are not needed in the following. 

To compare with the result of Sec. lIII A| we should now 
calculate, for example, the component gp]j'^^ ■ However, 
expecting the decay exponent to be independent of a and 



a 9pq 



(a, a) 



Due 



/?, we consider the simpler case Tr[gp^ 

to the orthonormality ^""^'^""''* ~ ^mm the depen- 
dence on the vector components then drops out. Thus 
for p 7^ q 



where we defined 



such that E^ 



.(3,5) 



(35) 



ikd 



[e(3,5)(fc)]2 ^ A±B(k). Defining now z = 
the integral can be turned into a contour integral 



around the contour Izl = 1 



2EA 



dz 



{z - z+){z - z^y 



(36) 



where the poles z± are determined from the equation 
^2 _ [4^2 _ (^2 _^ 3^2)2](8^^3)-i^ 1 = 0. They are 

given by 



z± 



AA^ ~ {if -f 37^)^ 
167773 



± 



4^2 - (772 + 3^2)2 



167773 



(37) 



such that = l/z_, and we choose the signs so that 
z_ is inside the contour \z\ = 1. In addition to this, 
assuming that p < q, there is a pole of order q — p at 
z = 0. The integral can then be evaluated using residue 
techniques, with the result 



2EA 



(38) 



This leads to an exponential decay of the propagator with 
growing q — p > 0, when E is off- resonant (in which 
case z± are real- valued) . Using this result, we can give 
an estimate for the decay of the transmission function 
[Eq. ^] through a finite chain of length N by replacing 
Gi,MN with Tr[g^j^]/M. This yields 



r(i?)~[z+(i?)]-2A^ = e-2^'"[^+(^)l, 



and thus the exponent 



P{E) = 2ln[z+{E)]. 



(39) 



(40) 



It can be checked that this result is, in fact, equal to the 
result [Eq. (|26p] obtained for the finite chain. 

It is thus seen explicitly that the decay coefficient of 
the off-resonant transmission does not in any way depend 
on the couphng of the molecule to the leads. It should be 
kept in mind, however, that the relative position of Ep 
within the HOMO-LUMO gap depends on the electrode- 
lead coupling and the charge transfer effects. This infor- 
mation is still needed for predicting the decay exponent 
l3{Ep) of the conductance. 

The analytical results presented in this and the previ- 
ous section can be used for understanding the behavior of 
the transmission function upon changes in the parame- 
ters. For example, it should be noted that when 77 is made 
smaller, the band gap around E fa becomes larger, and 
at the same time the decay exponent P{E) grows. In 
this way, the conductance of a molecular junction can 
be controlled, for example, by introducing side groups to 
control the tilt angles (p between the phenyl rings . '^^'^'^ 
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FIG. 4: (Color online) (a) The coordinates of the carbon 
atoms in the direction z along the molecular wire. The left 
electrode is at 2; = and the length of a phenyl-ring unit is d. 
(b) Relative variation of the onsite energies for two different 
voltage profiles, A and B. The profile function P{z) describes 
how the harmonic voltage V{t) = V -\- VacCOS.{ujt) is assumed 
to drop over the junction, the voltage at z being given by 
V{z,t) = V{t)P{z). 



where are the distances of the carbon atoms from 
the left metal surface, and V{t) = V + VacCos{Lot). In 
case A, P{z) — inside the junction, while in case B 
P{z) = {L- 2z)/(2L), where L = Nd + d/3 is the dis- 
tance between the two metal surfaces. 

The profile B is more complicated, because the voltage 
ramp breaks the homogeneity of the wire. In this case 
the current must be calculated with the method outlined 
in App. [B1 In the case of profile A, however, the I-V 
characteristics can be calculated based on the knowledge 
of the zero-bias transmission function in the absence of 
light, t{E). As discussed in App.|3 the current is given 

bji^iiS^ 



00 2 r 

E h(f)] J dET{E + lhco) 



2e 

X [fUE] - fR{E)] 



(41) 



IV. PHYSICAL OBSERVABLES AND 
NUMERICAL RESULTS 



The low-temperature zero-bias conductance then takes 
the particularly simple formii^ 



In this Section we present numerical results based on 
our model. Throughout, we employ the "wide-band" ap- 
proximation for the lead self-energies, such that ^^(i?) = 
— iri/2 and Y,if {E) = —iTji/2^ with energy-independent 
constants T^^n . Furthermore we only consider the sym- 
metric case Tl = = F. First we briefly describe 
how we generalize the theory, as presented above, to take 
into account static and time-dependent voltage profiles. 
Then we concentrate on near-equilibrium (or "linear- 
response") properties, using as examples the conduc- 
tance, the thermopower, and the conductance enhance- 
ment due to light with low intensity and frequency. In 
this case, knowledge of the zero-bias transmission func- 
tion calculated above is sufficient, and we can discuss the 
length dependence of the transport properties in a simple 
way. After that we consider the dc current in the pres- 
ence of an ac driving field of more general amplitude and 
frequency, first concentrating on the case of infinitesimal 
dc bias, and finally on the I-V characteristics. 



A. Voltage profiles 

When considering finite dc or ac biases within a non- 
selfconsistent TB model that cannot account for screen- 
ing effects, one of the obvious problems is how to choose 
the voltage profile. Throughout the discussion, we shall 
refer to two possible choices, as depicted in Fig. HI They 
are in some sense limiting cases, and the physically most 
reasonable choice should lie somewhere in between. Pro- 
file A assumes the external electric fields to be completely 
screened inside the molecule, such that the onsite energies 
are not modified, while B corresponds to the complete ab- 
sence of such screening. In both cases, we can write the 
time-dependent onsite energies as Cp 



Gdcia,Lu)^Go J2 ["^'(f)] ^(EF + inw). (42) 
I— —00 

Here I indexes the number of absorbed or emitted pho- 
tons, Ji{x) is a Bessel function of the first kind (of 
order I), and a = eVac/fi'-^ is the dimensionless pa- 
rameter describing the strength of the ac drive. Note 
that Gdcia,LU = 0) = Gdc(a = 0,uj) = Gqt{Ef) = 
Gdc- Equation (j4ip may equally well be written in the 
form^iS 



1 — — OG 



Io{V + 2lhLu/e), (43) 



where Io{V) is the I-V characteristic in the absence of 
light [Eq. Below, the results from these formulas are 
compared to the numerical results for profile B. 

In Fig. O we plot the zero-bias transmission functions 
for wires with TV between 1 and 7. Notice that the four 
energy bands numbered 3-6 in Eq. p2|) are all visible, 
being separated by the HOMO-LUMO gap at E/'f « 
and the additional gaps at E/^ « ±1.7. Here we use the 
parameters r/7 = 5.0, (p — 40° (i.e. 77/7 « 0.77), and set 
the Fermi energy to Ep/'y = —0.4. These values are close 
to those used in Ref . |3y, where they were extracted from a 
fit to results for gold-oligophenylene-gold contacts based 
on density- functional theory (DFT). We shall continue 
to use them everywhere below. A DFT calculation for 
the HOMO-LUMO splitting of benzene, together with 
the results preceding Eq. yields the hopping 7 « 3 
eV. The length of a phenyl-ring unit is approximately 
d — 0.44 nm, and the largest ac electric fields Vac/L 
considered will be on the order of 10^ V/m. The photon 
energies tku will mainly be kept below the energy of the 
HOMO-LUMO gap of the oligophenylene. 
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the zero-bias transmission function t{E) as 



28.41.58.59 



FIG. 5: (Color online) Transmission functions for the 
oligophenylene wires with lengths = 1, 3, 5, 7. The param- 
eters are T/y — 0.5, if = 40°, Ef/"/ ~ —0.4, as discussed in 
the text. 
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FIG. 6: (Color online) Dependence of observables on the num- 
ber of units N: (a) Conductance, (b) Seebeck coefficient, and 
(c) the light-induced relative conductance-enhancement. The 
circles correspond to values extracted from the r{E) function 
[Fig. [5], using Eqs. dMl, (gS]), and (ge]). The red lines cor- 
respond to the simple order-of magnitude estimates of Eqs. 
(|47)) . (|48)) . and ((49)) . with the analytically calculated P{E). In 
(c) the crosses (x for profile A and + for profile B) show nu- 
merical results with the finite values a = 0.5 and huj/'y = 0.05 
(see Sec UVCj) . 



B. Near-equilibrium properties 

Let us start by illustrating the usefulness of the ana- 
lytical results of Sec. IIIII with a few examples. We con- 
centrate on low temperatures and small deviations from 
equilibrium. In addition to the linear-response conduc- 
tance 



(44) 



we shall consider the thermopower, or Seebeck coefficient. 
At low enough temperature T, this is given in terms of 



S = 



r'{Ep) 



3e t{Ef) 



(45) 



where prime denotes a derivative. Thus it measures the 
logarithmic first derivative of the transmission function 
at E — Ep. The sign of this quantity carries informa- 
tion about the location of the Fermi energy within the 
HOMO-LUMO gap of molecular junction."*! The third 
quantity we shall consider is the photoconductance. In 
the limit a ^ 1 and hio/"f 1 we can expand t{E) and 
the Besscl functions in Eq. ([42]) (see App. [X| to lead- 
ing order in these small quantities, yielding Gdc('^) — 
GoT{EF)+Go{anuj)'^T"{EF)/l6. Defining then the light- 
induced conductance correction AGdc{^) = Gdd^) — 
Gdc{uj = 0), where Gdc{i^ = 0) = Gdc = Got{Ef), the 
relative correction becomes 



AGdcia,u) _ {anujfT"{EF) 



Gdc 



16 t{Ef) 



(46) 



We thus see that this quantity gives experimental ac- 
cess to the second derivative of the transmission function 
ai E = Ep- Note that in this approximation, which 
can be seen as an adiabatic or "classical" limit, the 
conductance correction depends only on the driving field 
through the ac amplitude Vac = afioj/e. 

As discussed above, it is reasonable to assume that for 
large enough N , the transmission function t{E) satisfies 
the exponential decay law 



t{E) ^ C(£;)e-''(^)^ 



(47) 



at the off-resonant energies E k, Ep- Let us furthermore 
assume that C{E) is only weakly iJ-dependent. Then it 
is clear that the Seebeck coefRcient will have the following 
simple linear dependence on N (Refs. | 28ll36 ' 



S cx t'{Ef)/t{Ef) - -I3'{Ef)N. 



(48) 



In contrast, the light-induced conductance correction sat- 
isfies a quadratic law 



AGdc{uj)/Gdc (X t"{Ef)/t{Ef) 



-I3"{Ef)N +[l3'{EF)fN^. 



(49) 



Deviations from these laws can follow from the energy- 
dependence of C{E). 

In Fig. [S] we demonstrate these length dependences 
within our model for the oligophenylene junctions. The 
circles connected by lines show the results based on the 
transmission functions of Fig. O using Eqs. (|44| . (I45p . 
and (Uni)- The separate solid lines are the estimates of 
Eqs. (gT]), (gS]), and (gH), based on the analytic result for 
/3{E). The result for AGdc{^^)/Gdc is furthermore com- 
pared with some example results for finite a and w, using 
a = 0.5 and huj/^ = 0.05 (see below). Although Eq. (H^ 
was derived above by assuming the profile A, the result 
appears to be rather well satisfied for profile B as well. 
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FIG. 7: (Color online) Zero-bias conductance for diflterent 
driving frequencies and driving strengths a = eVac/huj. 
Panels (a)-(d) are for A'' = 1, . . . , 4. The solid lines correspond 
to profile A, and the dashed lines to profile B. The lower pair 
of curves is for a = 0.5, and the upper pair for a = 2.0. 
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FIG. 8: (Color online) Dependence of the conductance on 
A^. Circles represent the conductance in the absence of light, 
while the squares are for light with Tiu)/^ = 0.5 and a — 1.5. 
The solid line is for profile A and the dashed line for profile 
B. 



for the two voltage profiles A and B. This is seen in Fig. 
[71 where we show Gdc{i^) for = 1, . . . , 4 as a function 
of uj for two values of a, and for both profiles. The re- 
sults for profile A again follow from Eq. (|^^ . but the 
results for B require a more demanding numerical calcu- 
lation (see App.[B]). In both cases the effect of light is to 
increase the conductance considerably. The physical rea- 
son is that the photoassisted processes, where electrons 
emit or absorb radiation quanta, brings the electrons to 
energies outside of the HOMO-LUMO gap, where the 
transmission probability is higher. This happens when 
hu; exceeds the energy difference between the Fermi en- 
ergy and the closest molecular orbital, in this case the 
HOMO. The main difference between the two profiles is 
that in case B, the sharp resonances at some frequencies 
are smeared out, and thus the light-induced conductance 
enhancement tends to be smaller. The increase can still 
be an order of magnitude or more. 

The dependence of this effect on the length of the 
molecule is still illustrated in Fig. [U where the conduc- 
tances in the absence of light and in the presence of light 
with hu)/"f = 0.5 and a = 1.5 are shown as a function 
of A^. While the conductance in the absence of light has 
a strong exponential decay, in the presence of light this 
decay is much slower. For profile A the conductance ac- 
tually oscillates periodically, while in the case of profile 
B the oscillations are superimposed on a background of 
slow exponential decay. In the DFT-based results^ the 
oscillations were not present, or at least not visible for 
the cases A^ — 1,...,4 considered there. Indeed, they 
are likely to be artifacts of the our TB model that ne- 
glects all other than 7r-orbital contributions, as well as 
uses the wide-band approximation. 

The results of Fig. [8] can also be stated in terms of the 
relative conductance-enhancement AGdc{'jj)/Gdc- For 
large a and w, the increase of this quantity with A^ is 
exponential for both profiles A and B. This should be 
contrasted with the quadratic behavior for small a and 
LO [Eq. (HI])]. Thus, the fact that the results indicated 
by the crosses in Fig. [6] exceed the result of Eq. (|46)) is 
understandable. 



D. Current- voltage characteristics 



C. Zero-bias conductance at finite drive 
frequencies and amplitudes 

Next we consider the zero-bias photoconductance 
Gdci'-^) for light whose frequencies and intensities are 
not restricted to the adiabatic limit. We have discussed 
this case previously, based on DFT results for gold- 
oligophenylene-gold contactsi^ There, however, the anal- 
ysis was based solely on the simple formula of Eq. (|42)) . 
Here we show that those results are not expected to 
change in an essential way within a more refined theory, 
since the results of our TB model are not very different 



Finally we discuss the effects of light at finite voltages 
V. Let us first consider the properties of the I-V charac- 
teristics in the absence of light. Examples are shown in 
Fig. EJa) for the case N = 5. They consist of consecutive 
steps, which appear every time a new molecular level 
comes into the bias window between fi^ and fiji. These 
steps are seen as the peaks 1 and 2 in the differential con- 
ductance dl/dV shown in Fig.[9^b). The first one occurs 
roughly at the voltage Vi = 2{Ep — Ehomo)/^, where 
Ehomo is the energy of the HOMO. The factor 2 arises 
from the symmetric division of the voltages with respect 
to the molecular energy levels. In the case of profile B, 
the currents tend to be smaller than for profile A, but 
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FIG. 9: (Color online) (a) I-V characteristics (A'' — 5) with 
and without light for profiles A (solid lines) and B (dashed 
lines). Results in the presence of light with a — 1.5 and 
Twj/'y = 0.075 are indicated with an arrow, (b) The corre- 
sponding differential conductances, (c) Same as (b) , but con- 
centrating on the low-bias regime and on a logarithmic scale. 
The vertical dashed lines indicate the approximate positions 
of the main peak and the light-induced side peaks. They are 
all separated by 2Tujj/e in voltage. 




FIG. 10: (Color online) The voltage-dependent transmission 
function at three voltages for the wire with N = 5 and profile 
B. For profile A the result is independent of voltage and equal 
to t{E,V = 0). 



the current steps occur at roughly the same voltages. It 
should also be noticed that for profile B, a small negative 
differential conductance is present following some of the 
steps. The origin of this is the localization of the molec- 
ular eigenstates due to the dc voltage ramp, which sup- 
presses the transmission resonances. This can be seen 
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FIG. 11: (Color online) Same as Fig. |9jc) but for wires with 



in the voltage-dependent transmission functions t{E, V) 
in Fig. [ini 

In the presence of light, the step structure of the I-V 
curves is modified. For profile A, the results follow simply 
from Eq. (|4T|) or (|43)) . but for profile B a fully numeri- 
cal treatment is again needed. In Fig. [9] the results for 
a = 1.5 and Tiuj/^ = 0.075 are shown as the curves indi- 
cated with arrows. In Fig.[5I^a) it is seen that the current 
for voltages below the steps is increased, and decreased 
above them. This removes the negative differential con- 
ductance present in the case of profile B. These changes 
are associated with the appearance of additional current 
steps. Here we concentrate only on the additional steps 
in the low-bias regime at voltages ^ ^ Vi , as the relative 
changes are largest there. Fig.[5]Jc) shows the differential 
conductance on a logarithmic scale in this voltage region. 
It can be seen that there are multiple extra peaks below 
the main peak, all of which are separated by voltages 
2fiu)/e from each other. These peaks are "images" of the 
main peak a,tV = Vi, and are easily understood based on 
Eq. . For profile B all the peaks are moved to slightly 
smaller voltages and their spacing is reduced, since finite 
voltages tend to also suppress the transmission gap (see 
again Fig. [10]). Notice that, in contrast to high dc bi- 
ases [Fig. |nUa,b)], in the low-bias regime [Fig. [5]^c)] the 
results depend only weakly on the choice of the voltage 
profile. Thus the predictions of the model appear to be 
robust. To observe the side steps, the radiation frequency 
should be large enough such that the steps are not "lost" 
under the broadening of the main steps. On the other 
hand, it should be small enough to have at least one step 
present. Thus, if the voltage broadening of the main 
step at = Vi is approximately Ai/e, then we require 
Ai <huj < Ep ~ Ehomo- 

Figure [11] additionally shows the low-bias differential 
conductances for TV = 1,...,4, with other parameters 
chosen as in Fig. [H^c) . It is seen that the effects of light 
quickly become weaker, as the length of the molecule 
decreases. In the case = 4, small side peaks are still 
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observed. Larger effects could be obtained by increasing 
the parameter a. 

Similar-looking additional steps are visible in the I- 
V characteristics of an extended-Hiickel model for xylyl- 
dithiol in Ref . i8|. Despite the differences in magnitudes of 
parameters, and slight asymmetries in the geometries, it 
is likely that some of those steps have essentially the same 
origin as explained above. However, the most striking re- 
sult in that reference was the overall order-of-magnitude 
increase in the current. 



V. CONCLUSIONS AND DISCUSSION 

In this paper we have studied a 7r-orbital tight- 
binding model to describe elastic and photoassisted 
transport through metal-molecule-metal contacts based 
on oligophenylenes. In contrast with simpler linear chain 
models that have previously been studied in great detail, 
our model describes a specific molecule, and its parame- 
ters can be directly associated with quantities obtainable 
from DFT simulations, for example. Models of this type 
can be of value in analyzing the results of more detailed 
ab-initio or DFT calculations,'^^ and in making at least 
qualitative predictions in situations where such calcula- 
tions would be prohibitively costly. 

We first showed that at zero voltage bias the model can 
be studied analytically in a similar fashion as the simpler 
linear chain models. In particular, we derived an expres- 
sion for the decay exponent of the off-resonant trans- 
mission function. We then discussed the length depen- 
dence of the dc conductance, the thermopower, and the 
relative light-induced conductance enhancement in the 
case of light with a low intensity (a) and low frequency 
(lo). The conductance enhancement was found to scale 
quadratically with length. For large a and u, the rel- 
ative enhancement increases exponentially with length. 
Finally it was shown, by numerical calculations, that the 
current-voltage characteristics are modified in the pres- 
ence of light by the appearance of side steps with a volt- 
age spacing 2hu!/e. We demonstrated that the predic- 
tions of the model are robust with respect to variations in 
the assumed voltage profiles. This provides further sup- 
port for our previous results on the photoconductanccr 

In our work, only symmetrical junctions with sym- 
metrical voltage profiles were studied. Asymmetries can 
modify our results through the introduction of rectifica- 
tion effects?^ and can change the positions of the light- 
induced current steps. The experimental observation of 
additional steps with a spacing related to the frequency 
of the light would nevertheless provide more compelling 
evidence for the presence of photoassisted transport than 
a conductance enhancement alone. The latter can also 
have other causes<^ 

We note that the light-induced current steps are similar 
to the steps observed in current-voltage characteristics 
of microwave-irradiated superconducting tunnel junc- 
tions, where they result from photoassisted quasiparticle 



tunnelingi^i^ In that case, the main difference is that 
the energy gap necessary for the effect is located in the 
macroscopic electrodes, while the transmission through 
the tunnel barrier depends only weakly on energy and 
voltage. As a result, the current steps have a voltage 
spacing of precisely fioj/e. These effects are exploited in 
the detection of microwaves in radioastronomy^^ Sim- 
ilarly, one may imagine properly engineered molecular 
contacts as detectors of light in the infrared or visible 
frequency range. 

In terms of our model, to increase the chances of ob- 
serving the light-induced current steps, the aim should 
be to minimize the broadening Ai/e of the first main 
current step at voltage Vi, and to maximize a. Also, a 
wire with a large enough Vi should be used. The broad- 
ening Ai is related to the sharpness of the transmission 
resonances, and thus to the length of the molecule and 
its couphng to the electrodes, described by F. A decrease 
of F, however, increases the importance of Coulomb cor- 
relations. Their effect on photoassisted transport has 
recently been discussed within simple modelsJ^i^ In- 
crease of a through the light intensity, in turn, increases 
the heating of the electrodes^ and the excitation of lo- 
cal molecular vibrations?^ These may affect the geome- 
try through thermal expansion'*^ and structural deforma- 
tions, but will also give rise to an incoherent component 
to the current!^ At high enough photon energies, also 
the direct excitation of electrons on the molecule may be- 
come important. The relaxation of such excitations due 
to various mechanisms (creation of electron-hole pairs in 
the electrodes, spontaneous light emission) should thus 
also be consideredi^ Also conformational changes of the 
molecule are possible.^ Finally, a proper treatment of 
screening effects on the molecule and in the electrodes, 
the excitation of plasmons, and their role in the field 
enhancement^ are other issues that should be studied in 
more detail. 

Of course, for the investigation of most of these is- 
sues, noninteracting models of the type presented above 
are not sufficient. Strong time-dependent electric fields 
may have effects that can only be captured by self- 
consistent theories taking properly into account the elec- 
tron correlations due to Coulomb interactions. These 
interactions may influence the electronic structure in a 
way that would, at least, require the parameters of our 
model to be readjusted in the presence of the fight. 
Even the geometry of the junction can become unsta- 
ble, and so it should in principle be optimized with the 
light-induced effects included. Time-dependent density- 
functional theory is showing some promise for the treat- 
ment of such problemSf^i^ In addition to DFT, also 
more advanced computational schemes are being devel- 
oped to handle correlation effects i^i^ A systematic in- 
vestigation of the optical response of metal-molecule- 
metal contacts, and thus the testing of the predictions 
of the simple models^i^i^iii^i^ remains an important goal 
for future research. 
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APPENDIX A: SIMPLIFIED FORMULA FOR 
THE TIME-AVERAGED CURRENT 

Consider the expression Eq. ([2]) for the time-averaged 

(or dc) current. The coefhcient t^j^I{E), for example, is 
the sum of the transmission probabilities of all transport 
channels taking the electron from energy E on the left 
to energy E + khuj on the right. That is, for A; > 
(k < 0) it describes electron transmission under the ab- 
sorption (emission) of k photons. Assuming the wide- 
band approximation and the voltage profile A, Eq. ^ 
can be written in the more transparent forms of Eqs. (|4ip 
and (|43|) . This can be demonstrated rigorously using the 
equations of App. [b1 but it is instructive to consider the 
following simpler derivation. The idea is the same as in 
the "independent channel approximation" of Ref. 0- 

For now, we allow the ac voltage drops at the L and 
R lead-molecule interfaces to be asymmetrical. Thus we 
define the quantities ul and an, satisfying a — ul — 
aR. Since for profile A there is no voltage drop on the 
molecule, electronic transitions only occur at the lead- 
molecule interfaces. Thus the transmission coefficients 
'^lili^) ^'^s giv<5n by 

OO 

E [Ji-k{aR)fT{E+lhcj)[Ji{aL)]\ (Al) 

/— — oo 

where [J;(q;l)]^ is the probability for absorbing (emit- 
ting) I photons on the left interface and [Ji_fe(afl)]^ the 
probability for emitting (absorbing) I — k photons on the 
right interface. The propagation between the interfaces 
occurs elastically at the intermediate energy E + Ihuj, 
according to the transmission function t{E). A similar 

expression holds for t]^'^(_E). Using these and the sum 
formula Y.'^^-oo[Jk{x)f = 1, Eq. Q leads to 

Op °° r 
I{V;a,oj) ^ — dET{E + lhLu) 

/— — OO 

X [[Ji {aL)f hiE) ~ [Ji {ai,)f fR{E)] . 

(A2) 

Equation (|4l]) follows by setting — a/2 and an — 
—a/2, and the equivalent form of Eq. follows by 

changing summation indices and integration variables. 
Similarly, other suggestive forms may be derived . '^'^^i^^ 



For a; ^ 1 and I > one may expand J±i{x) « 
{±x/2y/l\ - (±a;/2)'+V(? -|- 1)!. This can be used in 
the hmit a <^ 1, huj/-f <^ 1 discussed in the text. 



APPENDIX B: GREEN'S-FUNCTION METHOD 
FOR THE TIME-AVERAGED CURRENT 

Here we outline the Green-function method^i^ii^ used 
for obtaining the results for voltage profile B. Consider 
again the dc current of Eq. In the case of a har- 

monic driving field, it is reasonable to assume the exis- 
tence of time-reversal invariance, in which case we have 
the symmetry^ 

rl^kE)=rl,t\E + khco). (Bl) 

The current expression of Eq. (8) in Ref. 4 was derived 
under this assumption, and that result can be brought 
into the form of Eq. Using the notation of that 

reference,— the coefficients can be written 

T^^liE) = Tr^[GiE)rfiE)GHE)r^^\E)] 
ri'^iE) = Tr^[G{E)t['\E)GHE)tfiE)], 

where the hats denote the extended "harmonic" 
matrices^^ and Tr^^ a trace over them. In particular, G 
is the matrix for the retarded propagator 

GiE) = [{E -Hi)-W- tL{E) - tR{E)]-\ (B3) 

where H is the Hamiltonian of the wire in the ab- 
sence of voltage profiles. The matrix E is defined by 
[E]m,n = [E + mhuj)Sm,n^, where m and n are the har- 
monic indices. Using the wide-band approximation for 
the electrodes, the matrices T,x and are given by 

\ ' (B4) 
[ry]„i,Ti(-E) = Jm-i{o!x)Jn-i{ax)'i^x, 

with X = L,R and ul.b. — ±a/2. Here Sjf is the self- 
energy matrix of lead X (extended to the size of H), and 
Fx — — 2ImSx- The matrix W includes the effect of the 
profiles for the voltage V{t) = V+Vac cos(wi). If W{t) = 
Wdc + Wac cos{ujt) is a diagonal matrix consisting of the 
onsite energies e^\t), then 

[W] 

m.n — + 5m+l,n)- (B5) 

In this formalism, the time-reversal invariance amounts 
to G and f^'^'^ being symmetric, i.e. — A. 

Equation (|Bip can then be proved by using the 
relations [G\r,i+k,7i+k{E) = [G],n,n{E + khuj) and 
[r^x^]m+k,n+k{E) = [t^^'''^U^niE + khu). We notc that 
f is defined with a different sign of I than in Ref. 0. 
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